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Abstract 

An analytic calculation of the differential cross section for elastic Compton scat- 
tering on the deuteron at photon energies uj in the range of 25 — 50 MeV is pre- 
sented to next-to-next-to-leading order, i.e. to an accuracy of ~ 3%. The calculation 
is model-independent and performed in the low energy nuclear Effective Field The- 
ory without dynamical pions. The iso-scalar, scalar electric and magnetic nucleon 
polarisabilities ao and /?o enter as free parameters with a theoretical uncertainty of 
about 20%. Using data at w La b = 49 MeV we find a = 8.4 ± 3.0(exp) ± 1.7(theor), 
(3q = 8.9 ± 3.9(exp) ± 1.8(theor), each in units of 10 -4 fm 3 . With the experimen- 
tal constraint for the iso-scalar Baldin sum rule, ao = 7.2 ± 2.1(exp) ± 1.6(theor), 
(3q = 6.9=F 2.1(exp) =p 1.6(theor). A more accurate result can be achieved by: (i) better 
experimental data, and (ii) a higher order theoretical calculation including contribu- 
tions from a couple of so far undetermined four-nucleon-two-photon operators. 
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Albeit scalar electric and magnetic polarisabilities are fundamental properties of the neutron, 
there is a fair amount of uncertainty in their values for the lack of free neutron targets. Thus, 
one has to resort to indirect methods to extract neutron polarisabilities. For example, Rose 
et. al. depend on a particular treatment of the neutron as a quasi-free particle inside the 
deuteron and find for inelastic Compton scattering on the deuteron — > jnp thatQ er n ) = 
10.7^| 3 7 UJ. Since the strong electromagnetic field of a heavy nucleus could enhance small 
effects from polarisability contributions, neutron scattering off a heavy nucleus is another 
method employed. However, analyses of low energy neutron-lead scattering experiments give 
conflicting values: a {n) = 12.3±1.5±2.0, /5 (n) = 3.1^1.5^2.0 in Ref. @, and more recently as 
part of an investigation of additional nuclei = 0±5 in Ref. [|J. In contradistinction, the 
scalar proton polarisabilities are well-measured from Compton scattering off the proton [Q: 
«(p) = 12.1 ± 0.8 ± 0.5 and (3^ = 2.1 =p 0.8 q= 0.5, respectively. 

On the theoretical side, a self-consistent and controlled description of the polarisabilities 
is only encountered with the dawn of Chiral Perturbation Theory (xPT) || . Unfortunately, 
estimates of corrections beyond leading order (LO) in x?T are difficult due to unknown 
parameters which already enter at next-to-leading order (NLO) [§,0]. In this paper, we 
propose low energy Compton deuteron scattering as a probe for neutron polarisabilities. 
This is not a new idea and has fervently been advocated by Levchuk and L'vov, by Weyrauch, 
and by Arenhovel, Wilhelm and Wilbois, most recently in Refs. || |9], |TI]] ; for its history see 
Ref. [ll[ and references therein. That the ;\RT values for the iso-scalar polarisabilities are 



consistent with Compton scattering data on the deuteron was demonstrated in two Effective 
Field Theory (EFT) NLO calculations in extensions of to the few nucleon system, one 
using perturbative pions [12], one using non-perturbative pions |13| . 



The present calculation is unique in the sense that all the theoretical approximations 
are explicit and controlled, up to the order of the perturbative EFT calculation. The 
polarisabilities extracted from this calculation are independent of any modelling of the short- 
distance dynamics of nucleons inside a deuteron or any heavy nucleus. In particular, they 
are independent of assumptions about the pion-nucleon interaction, independent of the form 
of the inter-nucleon potential, manifestly independent of the photon coupling to nucleons 
and pions as extended objects and independent of the choice of regulator. The calculation 
is analytic and straightforward. To the accuracy claimed, the only experimental inputs 
are simple low energy observables: the binding energy and pole residue of the deuteron, 
the nucleon iso-vector magnetic moment and the nucleon-nucleon scattering length in the 
1 So channel. Our result will therefore - to the accuracy claimed - be reproduced by any 
"realistic" potential model, and can hence serve as cross check. However, at higher orders 
beyond those presented here, there are contributions from certain four-nucleon-two-photon 
operators representing short-distance physics that might not be reproduced accurately in a 
model calculation. In conclusion, the polarisabilities are obtained in a model independent 
way by a fit to data, and no explanation as to their origin is attempted. 

The letter is organised as follows: We start by outlining the framework of the calcu- 
lation and the power counting needed to determine the theoretical accuracy. Then, the 
amplitude and cross section is presented. Finally, various fits of the polarisabilities to data 



1 In this letter, we express polarisability values in units of 10 4 fm 3 
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are discussed. A more detailed presentation will be given in an upcoming article |L4j] . 



Our analysis is in the context of the few-nucleon EFT in which nucleons and photons are 



the only dynamical degrees of freedom [jig] . This formulation can be viewed as a systema- 
tisation of Effective Range Theory to include interactions with gauge fields, and relativistic 
and short-distance effects. Contributions from pions and other, heavier degrees of freedom 
are included at low energy through multiple-nucleon-photon contact interactions in pertur- 
bation. Since one cannot apply this EFT at energies where cuts in the amplitude are probed 
due to pion propagation, meson exchange etc., the natural breakdown scale A is the pion 



mass « 140 MeV. For details of the formalism, we refer to a recent review [16]. 

EFTs live from the fact that at low energies, a separation of scales exists: Like the binding 
momentum of the deuteron 7 = 45.7066 MeV, the photon momentum k and energy u — \k\ 
are much smaller than A. There is another "high-energy" scale, namely the iso-spin averaged 
nucleon mass M N = 938.92 MeV. Thus, each physical low energy observable at a typical 
momentum scale Q ~ 7 is expanded systematically in powers of two small parameters: ? 
and jjj-. For convenience, we formally identify A/Mn ~ Q/A, which is accidentally valid. 
Thus, relativistic corrections with a typical size r ) 2 /M% ~ (Q/M N ) 2 = (Q/A) 2 x (A/M N ) 2 
contribute at Q((Q/Af), i.e. N 4 LO, see Refs. 

The centre-of-mass photon energy u> and the conversion of relativistic photon energies 
to non-relativistic nucleon energies introduce new scales: uj and \JujM^. Therefore, it is 
convenient to separate the calculation into two different energy regimes |[1^]. In regime I, 
uj < 7 2 /Mat, and ^M n uj ~ 7 ~ Q. We concentrate on regime II, i.e. higher photon energies 
on j in which factors of ^uMjy ~ A are summed to all orders. In addition, we make 
some integral approximations that allow us to effortlessly move from one regime to another 
with one simple expression pT2[] . Details of the power counting, calculation and validity of 



the approximations will be provided in a future publication |14j . 

In the Weyl gauge A = 0, the amplitude for Compton scattering can be written as 



• e2 

iA = 



2M, 



N 
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representing the scalar S, vector V and tensor T coupling of a photon with incoming (out- 
going) polarisation e 7 (eM with a deuteron of incoming (outgoing) polarisation e* d (e^). 

|- = a = 1/137. The power counting reveals that the vector and tensor amplitudes start 
contributing to the unpolarised cross section only at N 4 LO. Hence the most dominant term 
is the scalar one. The centre-of-mass scalar amplitude S is conventionally written as 

S = F(u, 6) e 7 ■ e 1 ; + G(u, 6) {% x e y ) ■ (t' x ^) , (2) 

where F and G represent the scalar electric and magnetic amplitudes, and 9 is the angle 

between the directions of the incoming and outgoing photon momenta, k and k' . The 
familiar Thomson limit is reproduced by F{u = 0) = 1 and G(uj = 0) = 0. 

For F and G, an expansion in powers of Q/A exists, F = F (0) + F^ + ... etc., where 
the superscript denotes the contribution of order (Q/A) n . We keep all N 2 LO, i.e. (Q/A) 2 , 
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corrections compared to LO, 0((Q/ A) ). One such effect is the deuteron wave function 
renormalisation, introducing factors of Z^ — 1. = 1.690(3) is the residue of the scattering 
amplitude at the deuteron pole and corresponds to the normalisation of the deuteron wave 
function at large wave length. — 1 is treated as 0(Q/A), see Ref. for details. Higher 
order corrections such as relativistic effects and 3 Si- 3 Di mixing are irrelevant for the present 
N 2 LO calculation 
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|18|1 . Assuming that all contributions not taken into account are of 
natural size, we expect our result to be accurate to order (Q/A) 3 < 0.03. This error estimate 
lies at the heart of a model-independent, controlled approximation provided by EFT. 

Now for details of the calculation, performed in the centre-of-mass frame. The contribu- 
tion from the seagull diagram in Fig. [I] (a) scales as 0((Q/A) ), which constitutes a LO 
contribution in both energy regimes. We get 
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Figure 1: Contributions to Compton scattering on the deuteron (represented by the circles) 
up to order (Q/A) 2 . Solid lines: nucleons; wavy lines: photons; squares: four-nucleon 
contact interactions needed to reproduce the correct deuteron residue; large disk: interactions 
via nucleon polaris abilities; dots: magnetic (Fermi) interactions; dashed double line: 1 So di- 
baryon. Graphs with a permutation of external lines or vertices are not displayed. 

Figures ffl (b) and (c) contribute at NLO in both regimes, and up to N 2 LO 

F b+C = -{Z d - 1) , G b+C = . (4) 

Next, we consider the contributions from the iso-scalar, scalar electric and magnetic nucleon 
polarisabilities with Lagrangean £ po i = 2n(a E 2 + (3 B 2 )N^N, where a := (a^ + a^)/2 
and p := (p® + (3 {n) )/2. The polarisability dia grams in Fig. [1] (d) contribute at different 
orders in the different energy regimes: 
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The nucleon as point-like particle can only be polarised by its pion cloud, and at LO in 
XPT a ,(3 ~ a/(87rf 2 rn n ), i. e . numerically 2 (a , A)) M N /a ~ 1/A 2 with / ff = 131 MeV « 
m„- p. Thus, in regime I these diagrams are suppressed by (Q/A) 6 . In contradistinction, 
they contribute at N 2 LO, (Q/A) 2 , in regime II. 

The power counting for the diagrams in Fig. |l] (e) and (f ) is slightly more involved because 
of factors of \JioM^ appearing in the propagators of the loop integrals from the conversion of 
relativistic photon energies to non-relativistic nucleon energies in the intermediate nucleon 
propagators. Some diagrams contribute at different orders in the different regimes. For 
example, in regime II (cu ~ 7), Fig. [I] (e) is of order 7 /\JujM n ps Q/A and contributes at 
NLO to F. Its contribution to the magnetic part G is suppressed by (Q/A) 2 , i.e. N 2 LO. But 
for small photon energies u < r y 2 /M^, Siegert's theorem requires this diagram to contribute 
as strongly as the seagull, Fig. [I] (a), because of gauge invariance. Indeed, the contribution 
to F from Fig. [J] is 0((Q/A) ), i.e. LO in regime I; to G it starts at N 4 LO. We find for 
Fig. |l| (e) in regime II 
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where the approximations made for analytic results agree with the exact answer to within 
0.5%. 

The nuclear magnetic moment contributions from Fig. [I] (f ) are suppressed by factors of 
(kiQ/MnY in both the energy regimes. Here, the square of the iso-vector nuclear magnetic 
moment in nuclear magnetons, n\ = 5.536, is anomalously large, given an expansion pa- 
rameter of Q/A ~ 0.3. Numerically, Fig. [I] (f) is therefore expected to contribute at N 2 LO, 
and we choose to include it. The iso-scalar magnetic moment k ~ 0.4 does not enter at 
N 2 LO because it is of natural size. Fig. [I] (f3) contains intermediate two nucleon scattering 
through the 1 Sq channel with a scattering length of s °) = —23.714 fm ~ l/Q- Analytic 
and exact results for the momentum integrals agree again to better than 0.5%. 
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To summarise, in regime II (uj ~ 7) the electric and magnetic amplitudes to N 2 LO are 
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As noted above, the scalar polarisability contributions to Compton scattering are suppressed 
by (u/A) 2 . Up to N 2 LO, Compton scattering is therefore sensitive to the nucleon polaris- 
abilities (Xq, (3q only in regime II, where they enter as the only undetermined parameters, see 
(|5]) and (^]) . Thus photon energies uj ~ 7 ~ 45 Me V are most appropriate for extracting nu- 
cleon polarisabilities from a comparison of the differential cross section in the centre-of-mass 
frame 



da a 2 



2 (u + y/u 2 + 4M 2 )' 



|G| 2 ) (l + cos 2 #) +4Re[FG*]cos#] , (9) 



expanded to N 2 LO, with experimental data. Using the experimental values for the proton 
polarisabilities |4[], we can then determine the neutron polarisabilities. 

We now proceed to determine the polarisabilities from Compton scattering data at 
uj Lah = 49 MeV ||g. If the LO %PT result a = 12 A, p = ±a [g is taken at face 



value, one might assume that (3 « at N 2 LO. A x 2 fit gives ocq = 8.1 ±4.6, with sizable la 
errors^. Estimates of the NLO corrections in %PT seem to favour the electric and magnetic 
polarisabilities to have about the same size: Resonance saturation suggests «o = 12.0 ± 2.5 
and /3o = 5.7±5.1 0, and treating the A as dynamical degree of freedom leads to ao = 16.4, 
Po = 9.1 [0. In the next step, we therefore fit both ao and Po as free parameters, resulting 
in oq = 8.4 ± 3.0 ± 1.7 and (3q = 8.9 ± 3.9 ± 1.8, again with large error bars from a two pa- 
rameter fit to four data points. The first error comes from the fit, the second is the estimate 
of the theoretical uncertainties. However, the well measured Baldin sum rule constrains 



«o + Po — 14.05 ± 0.67 [^TJ which is also comfortably close to the sum of the polarisabilities 
from the naive two-parameter fit. Utilising the sum rule, we find «o = 7.2 ±2.1 and hence 
Po = 6.9 =F 2.1 in a one parameter fit, with the best \ 2 of all the three procedures. Strictly 
speaking, the Baldin sum rule gives the sum of the polarisabilities only at zero photon en- 
ergy, while the polarisabilities are extracted at uj = 47.7 MeV. However, we assume that 
corrections are small |22| . Figure compares the three fits to scattering data. 



The calculation presented here is formally accurate to 0((Q/ A) 3 ), i.e. to within ps 3%. 
However, numerically the result might be slightly better as some contributions with sizes 
comparable to N 3 LO are already included in and (0) and the polarisabilities are by 
far the biggest N 2 LO contributions. Figure |2| compares the strengths of the diagrams, 
confirming the power counting and demonstrating convergence. We can for example leave 
out the hi\ contribution Fig. [l] (f) which acts like a large N 3 LO correction, see Fig. [| This 
enhances the electric polarisability using the Baldin sum rule to ao = 8.1 ± 2.1 and reduces 
Po accordingly. On the other hand, one might expect that wave function renormalisation 
effects on the polarisability diagrams are large because Z<i — 1 = 0.69 ... is large for a 
NLO effect. Still, the fit using the Baldin sum rule then induces only small changes in 
the polarisabilities, «o = 7.3 ±2.1. We therefore estimate the theoretical uncertainty of 



our extraction at lu ps 40 MeV to be about 20% and refer to ||T4| for further discussion. 
The experimental uncertainty is at present higher than the theoretical error because of the 
comparatively high error bars in the only available low energy Compton scattering data. 



^Unless otherwise stated, the errors are only from the fit and do not include the theoretical uncertainties. 
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Including the theoretical uncertainties in our calculation (« 20%) and in the value of 
the Baldin sum rule, we find a = 7.2 ± 2.1 ± 1.4 ± 0.7, /3 = 6.9 q= 2.1 =F 1.4 =f 0.7, 
where the first error comes from the fit, the second one is the theoretical uncertainty of our 
N 2 LO calculation, and the third comes from applying the Baldin sum rule. For curiosity, 
we note that these values are also well consistent with the two Urbana data points at 
^Lab = 69 MeV |2(| which give using the Baldin sum rule «o = 5.4 ± 2.2, (3 = 8.6 =F 2.2. 
We finally extract for the neutron polarisabilities ~ (5; 2) and (3^ ~ (16; 12), where 
the first number in parentheses comes from the free fit, and the second from the fit using 
the Baldin sum rule. We quote no errors since they are large mainly due to the large 
experimental uncertainties. Higher values for seem slightly favoured Hli |. 
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Figure 2: Left: EFT result for the differential cross section fitted to data 0JjJ . Dashed: 
°-o — Po = 0; dot-dashed: (3q = 0; dotted: olq and /3q fitted; solid: «o and (3q constrained 
by the Baldin sum rule. Right: Sizes of the (angle averaged) contributions Fi and Gi, 
normalised to the LO contribution of the seagull term: Gi\/\F^\. The values assumed 
for the polarisabilities are 12 and 5. The grey bars are meant to guide the eye in determining 
the size of a contribution in the Q/A expansion. 

The gap between our extraction of the polarisabilities and the results is well known, 
It can be traced back to the fact that the Compton scattering data at large 
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see e.g. 

angles is systematically more and more enhanced over forward scattering as the photon 
energy is increased p3fl . This favours larger values of (3q. 

Although the scarcity of data at low energies is a big hindrance, the comparison of the 
present calculation to experiment shows that it is quite appropriate to determine neutron 
scalar polarisabilities from low energy Compton scattering. Our approach uses EFT and 
hence it is mo del- independent and allows one to estimate the theoretical accuracy as ~ 20% 
for the polarisabilities, assuming higher order contributions are of natural size. Future 
higher precision experiments e.g. at TUNL in the energy regime to ~ 25 — 50 MeV would 
be most useful. As one can see from Fig. |2j, at lower photon energies, the polarisability 
contributions become smaller than N 2 LO, while higher energies introduce large theoretical 
errors since corrections going like u/A are not suppressed sufficiently strong any more. The 
energy regime proposed is hence an interesting window of opportunity to determine nucleon 
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polar isabilities in a model-independent way without having to deal with pions as explicit 
degrees of freedom. Small angle scattering data constrains «o + flo and hence can serve as a 
cross check to the Baldin sum rule, whereas backscattering data probes ao — j3o and hence 
is needed for a better determination of (3q. In a future publication |TJ]], we will address 
details of the calculation as well as a more accurate estimate of N 3 LO effects. The only 
undetermined operators entering at N 3 LO are (C E E 2 + C M B 2 )(N T PiN)^ (N T PiN) with 
projections Pj onto the 3 Si channel. These operators contribute to residual electric and 
magnetic deuteron polarisabilities, i.e. to those not generated by polarisation of the two 
nucleons against each other or by polarising the nucleons themselves. 
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